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I.  INTRODUCTION 


For  many  years  maritime  cyclones  have  been  of  interest 
to  military  and  civilian  communities  froa  both  a  scientific 
and  a  practical  point  of  view.  Oceanic  lows  frequently 
develop  rapidly  and  produce  meteorological  phenomena  which 
can  adversely  impact  civilian  and  military  operations  in 
oceanic  and  coastal  regions.  Accurate  and  dependable  pre¬ 
diction  of  the  development  and  movement  of  oceanic  cyclones 
is  required  for  reliable  maritime  forecasts. 

Numerous  studies  of  middle-latitude  cyclones  have 
focused  on  continental  cyclones  (e.g.  Smith  (1973)  , 
Petterssen  and  Smebye  (1971),  Johnson  and  Downey  (1976), 
Chen  and  Bosart  (1977),  and  others).  This  is  primarily  due 
to  the  abundant  observational  data  which  is  readily  avail¬ 
able  from  continental  regions.  In  contrast,  relatively  few 
studies  have  investigated  maritime  cyclones  because  of  the 
limited  observational  data  from  oceanic  regions.  As  a 
result,  the  factors  influencing  maritima  cyclogenesis  are 
not  fully  understood. 

Several  studies  of  oceanic  cyclones  indicate  that  sig¬ 
nificant  differences  exist  in  the  physics  involved  in  conti¬ 
nental  and  oceanic  cyclogenesis.  Petterssen,  Bradbury  and 


=1 


rr 


Pedersen  (1962)  compared  the  development  of  cyclones  over 
the  western  North  Atlantic  to  that  of  cyclones  over  the 
eastern  United  States.  This  observational  study  indicates 
that  the  initial  development  of  oceanic  cyclones  appeared  to 
result  primarily  from  low-level  thermal  advection  while  con¬ 
tinental  cyclones  usually  developed  as  a  result  of  an 
upper- level  vorticity  maximum  approaching  a  low-level  baro- 
clinic  zone,  A  study  by  Nitta  and  Tamamoto  (1974)  also  sug¬ 
gests  that  maritime  cyclones  frequently  develop  in  the 
vicinity  of  an  ocean  baroclinic  zone  in  conjunction  with 
nearly  zonal  upper-level  flow.  Sanders  and  Gyakum  (1980) 
found  that  explosive  cyclogenesis  is  predominantly  a  mari¬ 
time,  cold-season  event  and  suggested  that  the  roles  played 

v 

by  cumulus  convection  and  air-sea  fluxes  in  cyclogenesis  are 
not  adequately  understood.  A  numerical  study  by  Sandgathe 
(1981)  indicates  that  air-sea  fluxes  play  a  significant  role 
in  the  evolution  of  oceanic  cyclones.  These  studies  and 
others  suggest  that  the  air-sea  fluxes  of  latent  and  sensi¬ 
ble  heat  play  a  significant  role  in  oceanic  storm 
development. 
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Current  operational  numerical  models  have  been  shown  to 
be  seriously  deficient  in  predicting  maritime  cyclone 
deepening  rates.  The  Sanders  and  Syakua  (1980)  study  indi¬ 
cated  that  the  National  H ete or o logical  Canter  (NHC)  primi¬ 
tive-equation  aodels  predicted  less  than  one  half  of  the 
deepening  rates  observed  in  explosive  oceanic  cyclogenesis. 
The  failure  of  numerical  aodels  to  predict  explosive  cyclo¬ 
genesis  is  not  entirely  the  result  of  the  Halted  observa¬ 
tional  data  froa  oceanic  regions.  Part  of  the  problea  lies 
in  modeling  vertical  exchanges  in  the  atmosphere.  The  upper 
troposphere  is  dominated  by  wave  structures*  while  vortex 
structures  dominate  in  the  lower  troposphere.  Host  cases  of 
upper  tropospheric  wave  amplification  are  accompanied  by  the 
development  of  anticyclones  and  cyclones  in  the  lower  tropo¬ 
sphere.  This  suggests  that  there  exists  a  strong  interac¬ 
tion  between  these  two  regimes.  The  smaller  size  of  oceanic 
cyclones  presents  another  problea.  The  space  scale  of  upper 
tropospheric  waves  is  typically  two  to  four  times  the  space 
scale  of  oceanic  cyclones.  &s  a  result,  it  appears  that 
numerical  models  adequately  simulate  the  large-scale 
processes  involved  in  synoptic-scale  systems,  but  do  not 
capture  the  smaller-scale  features  embedded  within 
synoptic-scale  weather  systems. 


Due  to  the  sparse  observational  network  in  oceanic 
regions,  few  quantitative  studies  of  maritime  cyclogenesis 
have  been  completed.  To  fill  this  void,  studies  of  numeri¬ 
cally  simulated  cyclone  development  are  needed.  Numerical 
models,  having  an  internally  consistent  "data"  base,  elimi¬ 
nate  errors  resulting  from  incomplete  or  inaccurate  observa¬ 
tions  from  data-sparse  regions  such  as  oceans.  In  cases 
where  there  is  a  favorable  comparison  between  numerically 
simulated  cyclogenesis  and  real  atmosphere  cyclogenesis,  the 
model  may  provide  a  greater  understanding  of  the  physics 
involved  in  cyclone  development  than  would  the  incomplete 
observations  of  the  real  atmosphere  storm.  K  disagreement 
between  the  model  and  the  real  atmosphere  would  not  neces¬ 
sarily  prove  to  be  of  little  or  no  value.  In  cases  where 
the  model  and  the  real  atmosphere  do  not  agree,  the  diagnos¬ 
tics  of  the  real  atmosphere  may  well  isolate  the  reasons  for 
the  disagreement  and  thus  lead  to  improvement  in  the  model. 

A  basic  feature  of  cyclogenesis  is  an  inward  mass  trans¬ 
port  (convergence)  in  the  lower  troposphere  and  and  outward 
mass  transport  (divergence)  in  the  upper  troposphere.  Asso¬ 
ciated  with  these  horizontal  mass  transports  is  an  upward 
vertical  mass  transport  within  the  cyclone.  Quantifying 


this  aass  circulation  can  provide  important  information  about 
the  large-scale  forcing  within  a  cyclone. 

Cyclones  are  vortices  embedded  in  a  planetary-scale 
pattern.  Consequently,  a  study  of  the  basic  properties  of 
cyclones  should  take  advantage  of  their  vortex  nature.  One 
such  approach  is  the  investigation  of  the  budgets  of  the 
fundamental  quantities  such  as  a  ass  and  circulation  within  a 
volume  moving  with  the  storm.  Quasi-Lagrangian  transport 
theory,  outlined  by  Johnson  and  Downey  (1975  a  and  b) ,  is 
particularly  suited  for  the  study  of  cyclones. 

In  this  thesis  the  mass  budgets  of  nuierically  simulated 
oceanic  cyclones  are  investigated  to  ascertain  the  basic 
aspects  of  cyclone  development  and  maintenance.  &  quanti¬ 
tative  study  of  the  aass  circulation  helps  describe  the 
intensity  and  physical  structure  of  the  cyclone  and  isolate 
vertical  variations  in  the  cyclone  circulation. 

This  study  follows  the  Ph.  D.  thesis  investigation  by 
Sandgathe  (1981).  The  budget  studies  were  applied  to  two 
storms  generated  by  a  version  of  the  OCLA  general  circula¬ 
tion  model,  which  includes  a  sophisticated  parameterization 
scheme  for  planetary  boundary  layer  (PBL)  and  free  atmo¬ 
sphere  interactions.  The  Sandgathe  study  demonstrated  the 


role  of  air-sea  flaxes  of  heat  and  moisture  in  initial 
cyclogenesis.  He  showed  that  cyclone  intensification 
resulted  when  sensible  heat  fluxes  were  resowed.  In  con¬ 
trast,  the  lack  of  noisture  fluxes  resulted  in  cyclone 
decay. 

In  this  thesis,  two  cyclones  were  studied — one  in  the 
southern  (spring)  heaisphere  and  one  in  the  northern  (fall) 
hemisphere.  Both  storms  are  characterized  by  rapid 
development  during  the  period  from  96  to  126  hours  of  simu¬ 
lated  cyclogenesis.  In  each  case,  the  central  surface  pres¬ 
sure  of  these  storms  decreases  at  a  rate  of  nearly  one  mb 
per  hour  during  this  30  hour  period.  Thu3,  these  storms  are 
good  examples  of  the  explosive  deepening  which  is  not 
adequately  handled  by  current  operational  dynamic  models. 

This  study  is  one  in  a  series  of  experiments  designed  to 
investigate  the  physical  processes  involved  in  extratropical 
cyclone  development.  The  primary  objectives  of  this  study 
are: 

•  to  study  the  synoptic  features  of  two  model-generated 
cyclones  and  compare  them  to  observed  cyclones; 

•  to  isolate  the  lateral  and  vertical  branches  of  the 

mass  circulation  during  different  stages  in  the 

evolution  of  a  model- genera  ted  maritime  cyclone;  and. 


•  to  compare  vertical  velocity  values  derived  through 
quasi-Lagrangian  diagnostics  with  the  vertical  velocity 
values  from  the  models  internally  consistent  ’•data" 
set . 
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The  atmospheric  model  used  for  this  study  is  briefly 
described  in  Chapter  II.  The  r a  suits  of  this  study  begin  in 
chapter  III  with  an  objective  review  of  the  synoptic  fea¬ 
tures  of  the  two  model  storms.  \  description  of  the  guasi- 
Lagrangian  diagnostics  technique  and  the  equations  used  in 
the  mass  and  circulation  budgets  are  presented  in  Chapter 
IT.  The  generalized  budget  equation  and  the  budget  volume 
coordinate  system  are  described  in  Appendix  B.  Chapter  7 
contains  the  analysis  of  the  mass  budget  results.  Conclu¬ 
sions  and  recommendations  for  further  studies  are  presented 
in  Chapters  71  and  711  respectively. 
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II.  A£flg££g££l£  32m  25S££IEII2S 


A.  MODEL  STRUCTURE 

The  sandgathe  experiments  used  the  Navy  Operational 
Ataospharic  Prediction  Systea  (NOGAPS)  aodel  which  is  being 
developed  by  the  Naval  Environaen tal  Pradiction  Research 
Facility  (NEPRF) .  The  NOGAPS  aodel  is  a  state-of-the-art 
version  of  the  UCLA  general  circulation  aodel.  A  brief 
overview  of  the  HOGAPS  aodel  is  presented  in  this  chapter* 
For  a  aore  detailed  description  of  the  NOGAPS  aodel  the 
reader  is  directed  to  Sandgatha  (1981).  The  UCLA  aodel  is 
described  in  detail  by  Arakawa  and  Laab  (1977). 

A  fine  resolution  grid  (2  degrees  latitude  by  2.5 
degrees  longitude)  within  a  60  degree  longitude  global 
sector  is  used  in  this  aodel.  This  grid  resolution,  which 
equates  to  a  222  ka  by  196  ka  grid  distance  at  45  degrees 
latitude,  is  sufficient  to  resolve  shorter  cyclone  waves. 

The  prognostic  variables  of  the  aodel  include  the  sur¬ 
face  pressure  ainus  50  ab,  tt  ,  the  horizontal  wind  velocity, 
V,  teaperature,  T,  and  specific  huaidity,  q.  Several  other 
prognostic  variables  (to  be  discussed  later)  are  associated 


with  the  PBL. 


The  aodel  uses  a  sigaa  vertical  coordinate  system 


V 


defined  as: 


where: 


P  *  pressure 


Ps  *  surface  pressure 

PT  a  pressure  at  the  top  of  the  aodel  ataosphere  (50  ab) . 


Thus  the  top  and  bottoay  Boundaries  of  the  aodel  are  a  »  0 
and  <T  *  1.  The  NOGAPS  aodel  consists  of  six  vertical  layers 
(Fig.  1)  .  In  this  depiction,  an  arbitrary  surface  pressure 
of  1000  ab  has  been  assuaed.  Variables  carried  at  the  layer 
aid- levels  are  T,  V,  and  g.  The  vertical  velocity,  b ,  is 
defined  at  the  layer  interfaces. 

The  aodel  eaployes  a  staggered  grid  (Arakawa  scheae  C) 
in  the  horizontal  (Fig.  2).  The  variables  defined  at  the 
center  points  are  T,  q,  and  na.  The  aaridional  wind  coapo- 
nent,  v,  is  carried  one  degree  north  and  south  of  the  center 
point  with  the  zonal  wind  coaponent,  u,  carried  1.25  degrees 
east  and  west  of  the  center  point.  For  diabatic  and  fric¬ 
tional  computations,  the  u  and  v  coaponents  are  averaged  to 
the  center  point.  A  saoothing  of  the  aaridional  coaponent 
adjacent  to  the  poles  is  inpleaented  in  theis  aodel. 

20 


PRESSURE  (MB) 


50- 


-  0 


0 


100 


200 


400 


600 


800 


jyjis 

ir<r 


\V,T,q 

Tfir 


VJL* 

ifa~ 


\VJ,q 


irar 


-J 

LU 

•3  > 
Ui 


1000 


.jyj.q  _ 

7 TO"  3  0 


6 


6 


Figure  1.  Vertical  distribution  of  sodel  large-scale 

prognostic  variables.  Pressure  values  of  signa 
levels  varv  with  the  surfacei pressure,  A  surface 
value  of  1000  nb  is  assuned  in  this  figure. 

(After  Sandgathe,  198  1) 
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LAYER 


The  area  used  for  this  study,  50  degrees  wide  extending 
from  pole  to  pole,  is  an  open-ocean  (i. e.  entirely  over 
water  with  no  sea  ice)  doaain  which  is  fixed  in  time  with 
two  zones  (centered  at  45N  and  35S)  of  maximum  north-south 
gradients  in  sea-surface  temperature.  Continuous  over-ocean 
propagation  within  the  model  is  guaranteed  by  imposing 
cyclic  boundary  conditions  in  the  east-west  direction. 
There  are  91  and  24  grid  points  in  the  north-south  and 
east-west  directions  respectively.  For  this  thesis,  zero 
degrees  longitude  has  arbitrarily  been  assigned  to  the 
western  boundary. 

The  PBL  paramererization  scheme  of  the  HOGAPS  model, 
which  follows  Deardorff  (1972)  and  Randall  (1976),  provides 
for  interactions  between  the  PBL  and  cumulus  and/or  stratus 
cloud  ensembles  at  each  grid  point.  Momentum,  heat  and 
moisture  fluxes  at  the  surface  are  diagnostically  determined 
using  the  values  of  T,  7,  and  q  from  the  dynamic  portion  of 
the  model  and  a  bulk  Richardson  number  which  is  based  no 
sea-surface  temperature.  The  cumulus  parameterization 
scheme  follows  Arakawa  and  Schubert  (1974)  and  Lord  (1978)  . 
Effects  of  water  vapor,  predicted  cloud  distributions ,  and 
large-scale  precipitation  processes  are  included  in  the 
radiational  heating  package  which  follows  Katayama  (1972). 
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Figure  2.  Horizontal  distribution  of  so$el  large-scale 
prognostic  variables.  This  distribution  fol 
Arafcawa  schese  C.  (After  Sandgathe,  1981) 


ollows 


B.  INITIAL  CONDITIONS 

The  sandgathe  study  utilized  initial  conditions  which 
are  representative  of  open-ocean  fall  and  spring  regiaes. 
The  initial  conditions  for  the  model  atmosphere  follow 
Sissons  and  Hoskins  (1977,  1978).  Hind  profiles  for 
the  fall  hemisphere  were  determined  using  a  polynomial 
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approximation  of  the  January  aean  zonal  wind  profile  at  30N 
froa  Oort  and  Basmussen  (1971).  The  faLl  hemisphere  winds 
are  salt  ip lied  by  a  factor  of  1.5  at  each  level  to  provide 
the  spring  hemisphere  wind  profiles. 

The  aodel  start  tiae  was  1  October.  Thus,  with  respect 
to  solar  fluxes,  the  northern  and  southern  hemispheres  are 
the  fall  and  spring  regimes  respectively.  Sea-surface 
temperatures  and  atmospheric  jets  are  also  representative  of 
these  seasons.  Zonal  jets  of  35  m  sec-1  and  52  a  sec-1, 
centered  at  45H  and  35S  respectively,  represent  possible 
conditions  for  these  regiaes.  is  shown  by  Simeons  and 
Hoskins  (1977,  1978),  the  existence  of  these  jets  guarantees 

the  growth  of  baroclinically  unstable  waves. 

% 

Sandgathe  used  the  conditions  listed  above  as  initial 
conditions  in  his  experiments  on  the  evolution  of  extratrop- 
ical  cyclones,  h  snail  amplitude  sinusoidal  perturbation  in 
wave  nuaber  six  was  added  to  the  meridional  wind  coaponent 
at  each  level.  The  maximum  value  of  the  perturbation  was  .5 
a  sec-1  at  45H/S.  The  perturbation  at  the  equator  and  the 
poles  was  zero.  The  associated  perturbation  in  the  surface 
pressure  field  was  approximately  one  millibar. 
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III.  £XNgm=  BIS£g£SI2£ 


A.  OTER  TIES 

An  initial  objective  of  this  thesis  is  a  coaparison  of 
the  synoptic  features  of  numerically  simulated  and  observed 
oceanic  cyclones.  This  discussion  mill  focus  on  the  signif¬ 
icant  features  in  the  model  "data”  fields.  Available  "data" 
consists  of  model  output  at  three-hour  intervals  from  hour 
96  to  hour  168  inclusive.  Except  for  the  1000  mb  heights 
and  surface  pressures,  all  fields  in  this  discussion  are 
generated  from  aid-layer  model  data.  With  the  exception  of 
the  wind  fields,  a  three  point  zonal  smoothing  technique  was 
applied  to  all  data  fields  represented  in  this  chapter.  The 
u  and  v  wind  components  were  averaged  to  the  height  grid 
points  for  the  wind  field  representations. 

The  synoptic  summary  of  the  two  model  cyclones  begins 
with  a  presentation  of  the  significant  features  in  the  lower 
troposphere  using  surface  pressure  and  temperature  fields, 
the  1000  mb  height  fields  and  the  1000  mb  to  mid-layer  4 
(approximately  1000-500  mb)  thickness  fields  (Figs.  4  to 
11)  .  This  is  followed  by  a  summary  of  mid-tropospheric 


features  (mid-layer  4,  approximately  530  mb)  of  the  height, 
temperature,  mind,  and  vorticity  fields  (Pigs.  12  to  19). 
Significant  features  in  the  aid-layer  3  (approximately  300 
mb)  height,  wind,  and  temperature  fields  are  described  in 
the  fourth  section  (Figs.  20  to  27).  This  chapter  concludes 
with  a  brief  comparison  of  the  development  of  the  model 
cyclcnes  and  observed  cyclones. 

Table  I  contains  a  list  of  synoptic  fields  relevent  to 
this  summary.  Fields  valid  at  96,  120,  144,  and  168  hours 
after  the  initial  time  are  presented.  With  respect  to  the 
synoptic  summary  and  the  figures  presented  in  this  thesis, 
the  indicated  times  represent  the  number  of  hours  of  numer¬ 
ical  simulation.  The  storm  tracks  of  the  two  numerically 
simulated  cyclones  used  for  the  budget  studies  of  this 
thesis  are  shown  in  Fig.  3. 

The  figures  presented  in  this  chapter  depict  two 
adjacent  60  degree  sectors  based  on  the  cyclic  continuity 
assumption  of  the  model.  Six  and  nine  degree  budget  volume 
circles,  centered  on  the  surface  pressure  minimum,  are 
depicted  on  each  chart.  Sadii  six  and  nine  were  chosen 
because  they  include  the  more  intense  portion  of  the 
cyclcnes  during  the  early  and  latter  stages  of  the  study 


inner  radius  and 


respectively.  In  following  discussions, 
outer  radius  refer  to  the  six  and  nine  degree  budget  radii 
respecti rely.  Both  the  northern  heaisphare  (NH)  and  south¬ 
ern  hemisphere  (SH)  cyclones  are  discussed,  relative  to  each 
figure,  to  point  out  differences  and  siailarities  in  their 
developaent. 


TABLE  I 

DATA  FIELDS  DEPICTED  IE  THE  SYNOPTIC  SUMMARY 


li?!i 

(OR  LAYER) 

lpfiSIIBISE 

LEVEL (S  ) 

FIELDS 

level  6 

surface 

pressure  (mb)  / 
temperature 

(1000  ab)/ 
aid-layer 

t 

4 

'  i88°ib 

1000  -  500  mb 
thickness/ 

1000  mb  height 

level  6 

surface 

pressure/ 

teaperature 

mid- layer 

4 

500  ab 

vorticity/ 

height 

aid- layer 

4 

500  ab 

temperature/ 

height 

aid-layer 

3 

300  ab 

tS’£g§!tuc,/ 

mid-layer 

3 

300  ab 

wind  vectors/ 
isotachs 
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Figure  3.  Store  tracks  of  two  aodel  cyclones.  Indicated 

locations  are  for  six  hour  intervals  froa  hour  96 
through  hour  168.  circle  radius  is  6°  latitude, 
a.  Northern  hemisphere,  b.  Southern  hemisphere. 


B.  LOWER-TROPOSPHERE  FEATURES 

Surface  pressure  and  temperature  fields  valid  at  hour  96 
are  shown  in  Fig.  4.  The  103  0  mb  height  and  1000  mb  to 
mid-layer  4  (approximately  1000  -  500  mb)  thickness  fields 
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rv 


for  the  saie  tiae  are  present9d  in  Fig.  5.  Rea*  thermal 

/ 

ridges,  at  32  (92)  and  57  (117)  degrees  longitude  respec¬ 
tively,  are  apparent  in  the  NH.  &  1011  tb  low  (cyclone  HI) 
is  located  in  the  theraal  ridge  near  32  (92)  degrees  longi¬ 
tude.  &  soaevhat  aore  intense  low  (cyclone  H2)  ,  with  a  cen¬ 
tral  surface  pressure  of  1007  ab,  is  in  the  theraal  ridg9 
near  57  (117)  degrees  longitude.  Initially,  cyclone  H2 

appears  to  be  the  aain  feature;  however,  within  24  hours, 
cyclone  N1  becomes  the  dominant  low  in  the  NH  sector.  k 
third  area  of  low  pressure  approximately  1000  na  northeast 
of  cyclone  HI  develops  very  little,  but  does  persist  through 
hour  156.  This  low  pressure  area  will  be  discussed  briefly 
later  in  this  section.  In  the  SH,  a  1003  ab  low  (cyclone  5) 
is  located  in  a  theraal  ridge  oriented  south-southwest  (SSW) 
to  north -northeast  (HHE)  and  centered  near  40  degrees  longi¬ 
tude.  Table  II  contains  the  central  surface  pressure  values 
of  the  three  cyclones  discussed  in  this  synoptic  summary. 
Cyclones  HI,  N2,  and  s  all  display  typical  lower- 
tropospheric  theraal  advection  patterns  with  warm  advection 
ahead  of,  and  cold  advection  behind  and  eguatorwardof  the 
developing  cyclones. 
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TABLE  II 


CEHTRAL  S0RPACE  PRESS  ORES  OP  THREE  MODEL  CYCLORES 

(in  ib) 


TIME 

CYCjj^RE 

CYCJjJME 

CYCLONE 

96 

1011 

1007 

1003 

99 

1008 

1002 

1001 

102 

1006 

1001 

998 

10  5 

1006 

1002 

996 

108 

1005 

1001 

994 

11 1 

1000 

998 

992 

114 

997 

996 

989 

117 

995 

998 

988 

120 

992 

998 

986 

123 

988 

997 

984 

126 

987 

997 

979 

129 

988 

998 

979 

132 

989 

1000 

979 

135 

989 

1001 

978 

13  8 

987 

1001 

976 

141 

989 

1004 

975 

144 

990 

1005 

975 

147 

988 

1006 

975 

150 

987 

1006 

976 

15  3 

988 

1008 

978 

156 

988 

1008 

980 

15  9 

985 

1006 

978 

16  2 

985 

1006 

979 

16  5 

985 

1006 

980 

16  8 

984 

1006 

981 
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P 
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* 


hours.  1000  ab  height  contours  (solid)  and 
thickness  contours  trashed)  are  m  aeters. 
a.  and  b.  as  in  Fig.  3. 
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Surface  pressure  and  temperature  fields  at  hour  120  are 
given  in  Fig.  6  with  height  and  thickness  fields  in  Fig.  7. 
During  the  period  from  hour  96  to  hour  123 ,  the  HH  lows  move 
approximately  700  nm  to  the  east  while  moving  only  slightly 
poleward.  Cyclone  HI  decreases  to  992  mb  and  becomes  the 
main  surface  low,  while  cyclone  N2  decreases  to  only  998  mb. 
Approximately  70%  of  the  surface  pressure  decrease  experi- 
enced  by  cyclone  HI  occurs  during  the  latter  half  of  this  24 
hour  period.  In  contrast,  more  than  half  of  the  surface 
pressure  decrease  experienced  by  cyclone  H2  occurs  during 
the  first  three  hours  of  this  period.  It  should  be  noted 
that  cyclone  HI  experiences  a  rapid  central  surface  pressure 
decrease  of  17  mb  between  hour  108  and  hour  123.  In  the  SH, 
cyclone  S  decreases  to  986  mb  at  a  nearly  constant  rate  of 
approximately  two  millibars  per  three  hours. 

The  thermal  ridges  associated  with  cyclones  HI  and  s  are 
more  pronounced  than  at  hour  96;  however,  thermal  advection 
at  the  surface  increases  very  little  during  the  period  from 
hour  96  to  hour  120.  By  the  end  of  this  period,  pockets  of 
war*  air  are  coincident  with  the  low  pressure  centers  in  the 
HH.  However,  cyclone  H2  does*  not  continue  to  develop  and 
thus  will  not  be  discussed  further. 
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Figure  6.  sane  as  Fig.  4  except  for  initial  tine  +  120 
nours. 
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Figure  7.  Saae^as  Fig.  5  except  for  initial  time  ♦  120 


The  analyses  for  hour  144  are  presented  in  Pigs.  8  and 
9.  During  the  period  froa  hour  120  to  hour  144 ,  cyclone  HI 
aoves  approziaately  700  na  to  the  east-northeast  (EHE)  while 
cyclone  S  aoves  approziaately  5 00  na  to  the  east-southeast 
(ESE).  Cyclone  HI  ezperiences  little  intensification  as  the 
central  surface  fluctuates  between  992  and  987  ab  with  the 
lowest  value  occurring  at  hours  126  and  138.  The  low  pres¬ 
sure  area,  initially  to  the  northeast  of  cyclone  HI,  inten¬ 
sifies  slightly  and  aoves  to  a  location  approziaatey  1000  na 
to  the  northwest  of  cyclone  HI  in  an  area  of  lower-tropos¬ 
pheric  wara  advection.  During  this  saae  period,  the  central 
surface  pressure  of  cyclone  s  decreases  to  975  ab  with  aore 
than  half  of  this  decrease  occurring  during  the  first  three 
hours  of  the  period.  kt  hour  1 32  (not  shown) ,  a  pocket  of 
wara  air  is  coincident  with  the  surface  center  of  cyclone  s. 
k  significant  feature  of  the  HH  thickness  field  is  the  north 
to  south  orientation  and  the  aaplitude  of  the  thermal  ridge. 
In  contrast,  the  SH  thermal  ridge  is  oriented  nearly  west  to 
east  and  is  eztreaely  elongated. 

Surface  pressure  and  temperature  fields  valid  at  hour 
168  are  presented  in  Pig.  10  with  the  height  and  thickness 
fields  given  in  Fig.  11.  During  the  period  from  hour  144  to 
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Figure  8.  Saie  as  Fig.  4  except  for  initial  time  ♦  144 
nours. 
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hour  156  hours,  cyclone  N1  reaains  relatively  steady-state 
with  a  central  surface  pressure  of  approximately  988  mb 
while  the  central  surface  pressure  of  cyclone  S  rises  5  mb 
to  980  mb.  Proa  hour  156  to  hour  168,  the  central  surface 
pressure  of  cyclone  SI  decreases  to  984  ab  while  cyclone  S 
reaains  nearly  constant  at  about  980  mb.  During  the  final 
24  hours  of  this  study,  the  SH  thermal  ridge  shifts  from  a 
north-south  orientation  to  a  northwest-southeast  orienta¬ 
tion.  However,  the  SH  thermal  ridge  maintains  its  west-east 
orientation.  The  1000  ab  height  contours  nearly  parallel 
the  thickness  contours  poleward  of  each  cyclone  so  there  is 
little  war a  advection  in  the  vicinity  of  these  cyclones; 
however,  cold  air  advection  continues  eguatorward  of  each 
cyclone.  During  the  period  from  hour  144  to  hour  162, 
cyclone  N1  moves  approximately  300  na  to  the  NNE.  It 
appears  that  cyclone  N1  moves  approximately  360  nm  during 
the  last  six  hours  of  this  study;  however,  this  apparent 
rapid  movement  may  actually  be  a  reformation  of  the  low  in 
the  vicinity  of  the  low  pressure  area  located  northwest  of 
cyclone  »1  at  hour  144.  Cyclone  S  also  exhibits  an  apparent 
rapid  movement  during  the  final  24  hours  of  simulation,  mov¬ 
ing  approximately  720  nm  to  the  ESE.  Approximately  500  nm 
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of  this  movement  occurs  during  the  period  from  hour  150  to 
hour  162  as  the  central  surface  pressure  increases  three 
millibars  to  979  ab.  As  in  the  NR,  it  is  possible  that  this 
apparent  rapid  movement  is  a  reformation  of  cyclone  S.  How¬ 
ever,  unlike  the  NH  case,  there  is  no  secondary  low  or  warm 
pocket  in  the  vicinity  of  the  new  location  of  cyclone  S. 

C.  HID- LE TEL  FEATURES 

In  this  section  mid-tropo spheric  features  are  discussed 
using  aid-layer  4  of  the  aodel  (approximately  500  ab)  . 
Hid- layer  4  height  and  teaperature  fields  valid  at  hour  96 
are  depicted  in  Fig.  12  with  the  vorticity  and  wind  fields 
given  in  Fig.  13.  A  weak  short  wave  within  nearly  zonal 
flow  is  apparent  in  each  hemisphere.  The  height  and  temper¬ 
ature  fields  are  nearly  in  phase  so  there  is  little  thermal 
advection  in  either  heaisphere.  Developing  cyclones  N1  and 
S  are  located  under  the  short  wave  ridges  in  the  respective 
hemispheres  with  a  weak  vorticity  maximum  located  300  to  400 
nm  to  the  west  of  each  cyclone.  Hi d-tropospheric  wind  flow 
patterns  indicate  weak  cyclonic  vorticity  advection  (CVA) 
over  each  cyclogenesis  area. 
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Figure  13. 


Mid-layer  4  absolute  .  vorticity  X?eP”lL  and  wind 
fields  valid  at  initial  tine  *  96  hours. 
Absolute  vorticity  (x  10— s)  contours  are  solid, 
a.  and  b.  as  in  Fig.  3. 
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Mid- layer  4  fields  valid  at  hoar  123  are  presented  in 
Figs.  14  and  15.  The  height  and  theraal  fields  reaain 
nearly  in  phase  with  weak  cold  advection  immediately  to  the 
west  of  the  center  of  cyclone  HI  and  warn  advection  approxi¬ 
mately  5  degrees  longitude  further  to  the  west.  During  the 
first  24  hours  of  this  study,  the  amplitude  of  the  NH  short 
wave  doubles  while  that  of  the  SH  short  wave  increases  by 
approximately  70%.  This  amplitude  intensification  is  nearly 
simultaneous  with  the  most  rapid  surface  development  which 
occurs  during  the  period  from  hour  96  to  hour  126. 

The  orientation  of  the  MB  short  wave  trough  is  northwest 
to  southeast  while  that  of  the  SH  short  wave  is  southwest  to 
northeast  as  each  cyclone  maintains  its  position  relative 
its  associated  short  wave.  The  cyclonic  vorticity  maximum 
in  the  NH  increases  in  magnitude  and  is  located  directly 
over  cyclone  HI  resulting  in  anticyclonic  vorticity  advec¬ 
tion  (&? A)  west  of  the  center  of  cyclone  HI.  The  increase 
in  the  vorticity  maximum  at  mid -layer  4  and  the  most  rapid 
surface  development  also  occur  simultaneously.  In  the  SH, 
the  cyclonic  vorticity  maximum  also  incraases  in  magnitude; 
however,  in  contrast  to  the  NH  case,  the  SH  vorticity  max¬ 
imum  does  not  move  over  the  storm  center.  Instead  the 
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vorticity  aaxiaua  aoves  to  the  east  as  cyclone  5  aoves 
southeasterly  to  a  position  poleward  of  the  vorticity  aax- 
iaua.  is  in  the  NH  case,  intensification  of  the  vorticity 
aaxiaua  occurs  siaultaneously  with  the  aost  rapid  surface 
develops ent.  Zero  vorticity  adv action  occurs  over  cyclone  s 
while  strong  CVA  occurs  to  the  north  and  northeast. 

Hid-tropospheric  fields  valid  at  hour  144  are  given  in 
Pigs.  16  and  17.  The  aaplitude  of  each  short  wave  continues 
to  increase  and  the  orientation  of  the  trough  axes  becoaes 
aore  west-east.  Extreaely  sheared  flow  is  evident  in  the  SH 
and,  to  a  lesser  degree,  in  the  HH.  The  short  waves  and 
their  associated  cyclones  aaintain  their  positions  relative 
to  each  other  during  the  period  between  hours  120  and  144  as 
the  height  and  theraal  fields  renain  nearly  in  phase. 
Although  aaintaining  its  position  over  cyclone  Ml,  the  vor¬ 
ticity  aaxiaua  in  the  SH  weakens,  having  reached  a  peak  at 
hour  132  (not  shown).  As  in  the  NH  case,  the  SH  vorticity 
aaxiaua  reaches  its  peak  at  hour  132  and  laintains  its  posi¬ 
tion  relative  to  cyclone  S.  In  each  case,  the  vorticity 
patterns  are  quite  elongated,  because  of  the  extreae  shear 
present  in  the  wind  field. 
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Figure  14.  Sue  as  Fig.  12  except  for  initial  time  ♦  120 
n  ours • 
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FJ  cure  15.  Same  as  Pig.  13  except  for  initial  time  ♦  120 
hours. 
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Figure  17.  Saee  as  Pig.  13  except  for  initial  tine  ♦  1«4 
hours. 


At  hour  168  (Pigs.  18  and  19) ,  cyclone  HI  is  positioned 
under  the  inflection  point  ahead  of  the  trough.  The  HH 
short  wave  continues  to  amplify  and  the  shearing  is  sore 
pronounced  as  the  axis  of  the  short  wave  is  nearly  east  to 
west.  Between  hours  144  and  168,  the  SH  upper  trough  fills 
leaving  a  split  flow  pattern  with  nearly  zonal  flow  to  the 
south  and  the  weakening  trough  to  the  north  of  46S,  the 
latitude  of  cyclone  S.  A  weak  closed  low  appears  approxi- 
aately  700  ns  to  the  northwest,  however  cyclone  S  fails  to 
develop  into  a  deep  cold  core  vortex.  The  vorticity  fields 
continue  to  weaken  and  become  more  elongated  during  the 
final  24  hours  of  this  study.  Broad,  relatively  unifors 
areas  of  cyclonic  vorticity  are  located  eguatorward  of  each 
cyclone.  These  vorticity  patterns  nearly  parallel  the  wind 
flow  indicating  little  vorticity  advection  in  either 
hemisphere. 

D.  UPPER-LEVEL  FEATURES 

Hid- layer  3  (approximately  3  00  mb)  height  and  tempera¬ 
ture  fields  valid  at  hour  96  are  given  in  Pig.  20  with  the 

wind  fields  given  in  Fig.  21.  Nearly  zonal  flow  is  evident 

> 

with  weak  cyclonic  flow  700  nm  to  the  west  of  each  cyclone 
and  weak  anticyclonic  flow  1000  nm  and  2300  nm  to  the  east 


of  cyclones  81  and  S  respectively.  The'  two  jets,  centered 
near  45N  and  35S,  are  nearly  zonal  and  sonewhat  stronger 
than  those  in  the  initial  field.  The  8H  and  SH  jet  maxima, 
approxiaately  45  ■  sec**1  and  55  m  sec-1  respectively,  are 
located  900  to  1500  ns  to  the  west  of  the  developing 
cyclones. 

Daring  the  first  24  hours  of  this  stady,  the  flow  imme- 
diately  to  the  west  of  each  cyclone  becomes  increasingly 
cyclonic  as  the  amplitude  of  the  trough/ridge  system 
increases.  it  hour  120  (Figs.  22  and  23)  the  cyclones  are 
positioned  slightly  east  and  poleward  of  the  base  of  the 
trough.  The  jet  maximum  in  the  MB  increases  only  slightly 
while  the  SH  jet  maximum  increases  to  over  60  m  sec-*. 
These  jets  are  located  approxiaately  600  nm  west  of  the  base 
of  the  respective  troughs. 

The  upper-level  charts  valid  at  144  hours  are  shown  in 
Figs.  24  and  25.  In  the  NH,  cyclone  81  is  located  ahead  of 
the  trough  under  nearly  meridional  flow  with  a  jet  maximum 
of  over  40  m  sec-1  located  approxiaately  5 00  na  to  the 
southeast  and  a  secondary  jet  maximum,  also  greater  than  40 
m  sec-1,  located  near  the  top  of  the  ridge,  approximately 
600  nm  to  the  north  of  cyclone  81.  In  the  SH,  there  are  two 
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nd  temperature  fields  yalid 
at  initial  time  ♦95.  Height  contours  (solid) 
are  in  decameters. ,  Isotherms  (dashed)  are  in 
°C.  a.  and  b.  as  in  Fig.  3. 
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jet  Bazina  (each  greater  than  SO  a  sec-1)  associated  with 
cyclone  S--one  behind  the  trough  and  one  in  the  base  of  the 
trough.  In  each  heaisphere  the  flow  ahead  of  the  trough  is 
basically  aeridional  and  the  short  wares  are  continuing  to 
increase  in  aaplitude. 

The  aid-layer  3  wind  field  ralid  at  hour  156  (not  shown) 
is  characterized  by  several  saall  disjointed  jet  aaxiaa 
greater  than  40  a  sec**1  in  the  HH.  There  appears  to  be  no 
organization  to  these  areas  as  soae  are  ahead  and  soae 
behind  the  trough.  By  this  tine,  the  SB  trough  has  tilted 
back  on  itself  and  is  oriented  nearly  we3t  to  east  and  the 
flow  south  of  the  trough  is  now  zonal,  k  developing  cut-off 
low  is  indicated  by  the  extra ne  shear  in  the  SB  height 
field.  Siailar  directional  shear  is  present  in  the  SB,  alt¬ 
hough  not  to  the  extent  of  that  in  the  SB. 

Figs.  26  and  27  depict  the  aid- layer  3  patterns  at  hour 
168.  In  the  SH,  the  position  of  cyclone  SI  relative  to  the 
trough  renains  unchanged  during  the  final  24  hours  of  this 
study.  Bowever,  an  organized  jet  aaxiaua  of  over  40  m  sec**1 
to  the  west  is  approaching  the  base  of  the  trough.  This 
aaxiaua  is  a  result  of  the  cyclic  continuity  assuaption 
eaployed  by  the  aodel  and  is  the  sane  aaxiaua  which  was  in 
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the  downstream  ridge  (north  of  the  cyclone)  at  approximately 
hour  144.  The  flow  directly  oxer  and  north  of  the  cyclone 
HI  reaains  basically  aeridional  bat  shifts  somewhat  to  the 
northwest  producing  stronger  directional  shear  to  the  north¬ 
west  of  cyclone  HI.  In  the  SH,  a  jet  maximum  greater  than 
60  m  sec-1  is  centered  in  the  base  of  the  trough  750  nm 
eguatorw'ard  of  cyclone  S ,  which  is  located  under  a  region  of 
very  weak  flow  at  mid-layer  3.  A  weak  closed  low  at  mid- 
layer  3  is  located  approximately  700  nm  northwest  of  cyclone 
S,  almost  directly  oxer  the  weak  closed  low  at  mid-layer  4. 
The  flow  equatorward  of  each  cyclone  continues  to  be  zonal. 

E.  COUP  AH  ISO  H  TO  OBSEHVED  CrCtOHES 

i  superficial  analysis  of  model  cyclogenesis  might  gixe 
the  impression  that  these  storms  are  realistic  in  exery 
aspect.  &  closer  examination  rexeals  them  to  haxe  a  number 
of  unusual  features.  In  this  section,  the  realistic  and 
anomalous  features  will  be  summarized. 

Although  cyclone  S  is  more  intense,  the  overall 
development  of  both  model  cyclones  resembles  that  of  storms 
observed  in  the  Horth  Atlantic  by  Petterssen,  et  al.  (1962). 
The  model  cyclongenesis  begins  in  the  lower  troposphere  in 
the  vicinity  of  a  baroclinic  zone,  under  nearly  zonal  flow 
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Figure  27 


aloft,  and  with  no  significant  vorticity  advection.  The 
upper-level  cold  trough  and  the  low-leval  cyclone  develop 
simultaneously  while  aaintaioing  their  position  relative  to 
each  other.  In  addition,  the  prisacy  contribution  to 
developaent  is  provided  by  lower-tropospheric  thersal  advec- 
tion  (at  least  during  the  early  stages  of  development)  and 
not  upper- level  vorticity  advection. 

During  the  early  stages  of  cyclogenesis,  the  lower- 
tropospheric  features  are  very  realistic.  Near  the  surface, 
the  simulated  cyclones  resemble  the  classical  cyclone  model 
with  the  cold  air  on  the  west  side  of  the  cyclone  and  warm 
air  on  the  east  side.  Cyclones  Ml  and  s  are  nearly  symmet¬ 
rical  for  the  first  24  to  48  hours  of  this  study.  However, 
the  subsequent  developaent  results  in  a  very  elongated 
structure,  lit  hough  asymmetrical  cyclones  are  not  uncommon, 
the  extreme  elongation  of  these  cyclones  does  appear  to  be 
atypical. 

The  lower-tropospheric  thermal  patterns  are  also  real¬ 
istic  during  the  first  half  of  the  period  studied.  The  warm 
advection  in  the  vicinity  of  the  storm  during  the  early 
stages  of  developaent  and  the  warm  air  coincident  with  the 
storm  center  during  the  latter  stages  of  simulation  would 


tend  to  support  cyclogenesis.  However ,  the  elongated 
thickness  patterns  and  the  orientation  of  the  thermal  ridges 
in  each  heaisphere  daring  the  later  stages  are  unusual. 

The  Hitta  and  Taaaaoto  (1974)  observational  study  sup¬ 
ports  the  findings  of  Petterssen,  et  al.  (1962).  In  addi¬ 
tion,  they  noted  that  oceanic  cyclones  generally  predominate 
in  the  lover  troposphere.  The  development  of  the  two 
cyclones  used  in  this  study  is  generally  limited  to  the 
lover  troposphere  although  cyclone  S  developes  weak  closed 
lows  at  aid-  and  upper-tropospha  ric  levels.  These  features, 
however,  are  extremely  weak  and  do  not  appear  until 
hour  168.  It  is  not  surprising  that  cyclone  s  developes 
to  the  upper  troposphere,  however  the  vertical  coupling 
of  this  development  is  somewhat  unusual.  This  feature  will 
be  discussed  in  more  detail  latar  in  this  section. 

The  above  characteristics  of  the  model  storms  which 
appear  to  be  realistic  are  basically  associated  with  the 
lower  troposphere.  On  the  other  hand,  the  anomalous  fea¬ 
tures  ace  generally  associated  with  the  vertical  configur¬ 
ation  of  the  cyclones  and  the  structure  of  the  mid-  and 
upper-tropospheric  levels  of  the  model,  mid-layer  4  and 
mid-layer  3  respectively.  The  development  of  the 
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waves  in  the  aiddle  and  upper  troposphere  appears  to  be 
atypical.  In  particular,  the  orientation  of  these  waves  and 
the  resulting  shear  are  quite  unusual.  is  noted  earlier, 
troughs  in  each  hemisphere  eventually  tilt  back  on 
themselves  creating  extreme  directional  shear. 

The  height  and  thermal  fields  at  mid-layers  4  and  3  are 
nearly  in  phase  throughout  the  72  hours  of  this  study.  Fur¬ 
thermore,  a  close  examination  of  the  height  fields  at  these 
two  levels  reveals  them  to  be  nearly  identical  in  shape  and 
vertically  stacked  throughout  the  entire  period  of  this 
srudy.  This  is  a  reflection  of  a  self-amplification  process 
and  is  likely  related  to  the  lack  of  occlusion  and  shallow 
structure  of  the  mature  cyclones.  In  addition,  the  upper 
troughs  and  the  cyclones  maintain  their  positions  relative 
to  each  other  throughout  the  period  of  this  study.  These 
two  features  are  consistent  with  the  observations  of 
Petterssen,  et  al.  (1962);  however,  the  distance  between 
cyclone  S  and  the  weak  upper  lows  which  develop  in  the  SH 
after  168  hours  of  simulation  is  unusual.  The  closed  upper 
lows  are  vertically  aligned  but  they  are  nearly  700  nm  to 
the  northwest  of  the  surface  center.  By  the  time  the  upper 
lows  develop,  cyclone  S  has  reached  its  peak  intensity  (975 


mb  between  hours  141  and  147)  at.  the  surface  and  has  begun 
to  weaken  somewhat  (to  981  mb) .  One  might  expect  the 
upper-level  lows  to  be  more  intense  and  vertically  aligned 
with  the  surface  low  during  the  nature  stage  of  the  cyclone. 

The  development  of  the  patterns  in  the  mid-layer  4  vor- 
ticity  fields  also  appears  to  be  anomalous.  In  particular, 
the  location  of  the  vorticity  maxima  with  respect  to  the 
cyclones  (over  cyclone  N1  and  equatorward  of  cyclone  S) 
during  the  final  48  hours  of  this  study  is  unusual.  The 
extremely  elongated  regions  of  relatively  uniforn  vorticity 
which  parallel  the  wind  flow  also  appear  to  be  unusual. 

In  conclusion,  these  model  storms  are  characterized  by 
both  realistic  and  unusual  features.  The  realistic  features 
are  generally  in  the  lower  troposphere  with  the  atypical 
features  generally  in  the  middle  and  upper  troposphere.  The 
vertical  resolution  of  the  model's  upper  troposphere  could 
possibly  be  a  factor  in  the  development  of  these  unusual 
features,  although  all  aspects  of  these  storms  are  not  typ¬ 
ical  of  oceanic  cyclones,  they  are  realistic  enough  to  be 
useful  in  determining  the  mass  and  circulation  budgets  of 
numerically  simulated  cyclones.  Quasi-Lagrangian  diagnostic 
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techniques  are  described  in  the  next  chapter.  in  analysis 
of  the  mass  budgets  of  these  stores  is  presented  in 
Chapter  7. 


IV.  MtelzkhS&MSlkS.  2IAGN051IC  Iisggiflag 


A.  QOASI-LAGRANGIAN  DIAGNOSTICS 

A  primary  goal  of  this  thesis  is  to  describe  cyclone 
development  quantitatively  through  an  analysis  of  oass  bud¬ 
gets.  The  quasi-Lagrangian  diagnostics  approach  provides  a 
aeans  of  achieving  this  goal.  Nash  (1978)  deaonstrated  the 
validity  of  the  quasi-Lagrangian  diagnostics  approach  in 
budget  studies  of  observed  and  numerically-predicted 
cyclones.  Basic  quasi-Lagrangian  diagnostics  techniques  are 
presented  in  this  section  with  the  aass  budget  equation  pre¬ 
sented  in  the  following  section.  A  description  of  the 
generalized  quasi-Lagrangian  diagnostics  budget  equation  and 
the  budget  voluae  coordinate  systea  (Nash,  1978)  are 
presented  in  Appendix  B. 

Traditionally,  budget  studies  have  employed  a  stationary 
budget  voluae  which  is  positioned  so  that  the  feature  being 
studied  moved  through  or  remained  within  the  budget  volume 
during  the  period  of  the  study  (e.g.  Saith  (1973),  Kung  and 
Baker  (1975),  and  others).  Since  the  horizontal  extent  of  a 
stationary  budget  voluae  is  greater  than  that  of  a  cyclone 


scale  feature,  it  is  not  possible  to  isolate  a  cyclone  from 
the  environment  in  which  it  is  sowing.  In  addition,  budget 
studies  eaploying  stationary  budget  volumes  provide  a  quali¬ 
tative  estimate  of  the  large  scale  forcing  processes 
involved  in  cyclone  development  and  decay.  Therefore,  the 
use  of  a  stationary  budget  volume  is  not  particularly  suited 
to  budget  studies  of  cyclones. 

The  guasi-Lagrangian  technigue  is  ideal  for  budget  stud¬ 
ies  of  cyclones  because  the  budget  volume  moves  with  the 
velocity  and  direction  of  the  cyclone.  &  translating  budget 
volume  affords  a  method  of  distinguishing  between  the 
effects  of  storm  translation  and  actual  cyclone  development. 
Thus  a  translating  budget  volume  is,  in  fact,  an  instrument 
for  examining  the  fluxes  across  the  budget  volume  boundaries 
by  eliminating  the  distort  ion  that  results  from  the  movement 
of  the  storm. 

Budget  studies  which  employ  a  stora-f ollowing  budget 
volume  provide  a  view  of  the  interaction  between  the  smaller 
scale  cyclone  and  the  planetary  scale  environment  in  which 
it  is  embedded.  Through  the  transport  relation  and  other 
governing  equations,  the  quasi -Lagrangian  budget  approach 
provides  a  method  of  calculating  the  underlying  scale 


interaction  processes.  This  approach  emphasizes  the 
transports  of  sass  through  the  boundaries  and  the  changes  in 
the  circulation  within  the  budget  volume.  Quasi- Lagrangian 
diagnostics  incorporates  transports  into  and  out  of  the 
budget  voluae,  as  well  as  internal  sources  and  sinks  of 
these  basic  properties. 

The  quasi-Lagrangian  diagnostics  approach  has  been  used 
in  a  variety  of  cyclone  studies.  Table  II,  partially 
extracted  froa  Hash  (1978),  is  an  updated  list  of  several 
recent  budget  studies  which  applied  quasi-Lagrangian 
diagnostics  to  observed  and/or  simulated  cyclones. 


TABLE  III 

QUASI-LA  GRAM  61  All  DIAGNOSTICS  BUDGET  STUDIES 
(partially  extracted  froa  Hash,  1978) 


PBOPBHTT 

Available  Potential  Energy 

Circulation 

Hass 

Absolute  Angular  Hoaentun 
Kinetic  Energy 
Absolute  Angular  Hoaentua 
Hass 

Circulation/ 

Angular  Hoaentua 


RESEARCHES  (S) 

Spaete  1974 
Hash  1975 

Johnson  and  Downey  1976 
Johnson  and  Downey  1976 
Chen  and  Bosart  1977 
Hash  1978 
Roaan  1981 

Conant  1982 


As  mentioned  earlier,  quasi -Lagrangiin  transport  theory 
provides  a  diagnostic  framework  which  is  ideally  suited  to 
the  study  of  cyclone-scale  circulations.  In  addition,  the 
internally  consistent  "data"  base  of  numerically  simulated 
cyclones  limits  computational  arrors  to  those  resulting  from 
interpolation  from  a  rectangular  grid  to  the  circular  grid 
of  the  budget  volume  and  to  round-off  errors.  With  the 
quasi-Lagrangian  diagnostics  approach,  a  cylindrical  budget 
volume  is  placed  around  the  model  storm  and  is  moved  with 
the  storm  center.  ?or  this  study,  the  surface  pressure 
minimum  is  used  to  indicate  the  storm  canter.  The  budget 
volume  extends  from  the  earth*s  surface  to  the  top  of  the 

atmosphere  with  the  axis  of  the  budget  volume  normal  to  the 

« 

earth's  surface.  The  budget  volume  is  actually  a  portion  of 
a  cone  with  its  vertex  at  the  earth's  center.  However, 
since  the  tropospheric  height  is  much  less  than  the  radius 

s 

of  the  earth,  the  budget  volume  radius  is  assumed  to  be 
constant  in  the  vertical.  The  budget  volume  for  this  study 
consists  of  the  seven  model  sigma  levels  as  defined  in 
Chapter  II.  The  radius  of  the  budget  volume  is  defined  in 
degrees  latitude  relative  to  the  storm  center.  Having 
defined  the  location  and  radius  of  the  budget  volume. 
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the  budget  of  any  basic  store  properties  (e.g.  aass, 
circulation,  angular  aoeentuu,  kinetic  energy,  and  eoisture) 
can  be  studied. 

B.  BUDGET  EQUATIONS 

For  this  study  the  guasi-La grangian  diagnostics  techni¬ 
que  is  applied  to  eass  budgets  of  two  numerically  sieulated 
cyclones.  The  eodel  experiments  defined  vertical  notion  to 
be  zero  at  the  upper  and  lower  boundaries  of  the  atnosphere. 
Thus  vertical  transports  across  these  boundaries  are 
eliminated  in  the  nass  budget  equations. 

Table  IV  contains  the  nass  budget  equation  in  the  signa 
coordinates  of  the  sodel.  In  its  nost  basic  fora,  the  gen¬ 
eralized  budget  equation  consists  of  three  terns:  lateral 
transport,  vertical  transport,  and  source  (or  sink). 

Since  there  are  no  sources  or  sinks  of  nass,  the  source 
(sink)  tern  of  the  generalized  equation  is  zero.  The 

renaining  terns  in  the  nass  budget  equation  represent  the 
lateral  aass  transport  (convergence/divergence)  and  vertical 
transport  (residual)  respectively.  When  the  limits  of  inte¬ 
gration  are  <r  *  0  and  a  *  1  (where  a  *  0)  ,  the  budget  resid¬ 
ual  tern  nust  be  zero.  If  not,  the  lateral  mass  transport 
tern  for  the  various  model  layers  are  adjusted  so  that  the 


integrated  residual  vanishes.  Lateral  aass  transport  and 
aass  tendency  are  coaputed  within  the  aass  budget  and  resi¬ 
dual  iabalances  are  then  coaputad  froa  the  tine  averages  of 
these  quantities.  Tertical  integration  of  the  layer  resid¬ 
uals  produces  a  vertical  aass  transport  which  is  then 
expressed  as  a  vertical  velocity  field. 


V.  BUDGET  A KALISIS 


A.  IHTHODUCTIOH 

The  analysis  of  the  quasi-Lagrangian  lass  budget  begins 
with  a  presentation  of  the  thereal  structure  of  the  cyclones 
in  teres  of  potential  teeperature  (Section  B)  .  This  is 
followed  by  a  discussion  of  the  pressure  (eass)  tendency 
(Section  C)  and  the  total  lateral  eass  transport 
(Section  D).  Vertical  eass  transports  are  discussed 

relative  to  the  nodel  vertical  velocities  (Section  E) .  In 
an  effort  to  determine  the  accuracy  of  the  vertical  velocit¬ 
ies  computed  in  the  mass  budget ,  budget  vertical  velocities 
are  compared  to  the  model  vertical  velocities  (Section  P) . 

Host  figures  presented  in  this  chapter  are  layer /time 
period  cross-sections.  Budget  calculation  values  at  levels 
and  "data”  times  mere  averaged  to  prodace  layer  and  time 
period  midpoint  values  respectively  for  the  time  sections. 
Exceptions  are  indicated  where  appropriate.  Inward  and 
upward  transports  are  represented  by  positive  values  in  all 


these  figures 


B.  THERMAL  STRUCTURE 

Area-averaged  aid-layer  and  "data"  tiae  potential  temp¬ 
eratures  are  presented  in  Pigs.  28  and  2  9.  The  vertical 
structures  are  siailar  in  the  two  cyclones.  The  patterns 
for  radii  six  and  nine  are  nearly  identical.  In  the  upper 
troposphere,  the  potential  teaperature  increases  slightly, 
while  decreases  are  seen  in  the  Biddle  and  lower  troposphere 
during  the  last  half  of  the  study.  During  the  period  of 
aost  intense  cyclogenesis  (hours  96  to  125)  only  saall  fluc¬ 
tuations  are  seen  in  the  area-averaged  potential  tempera¬ 
ture.  The  cooling,  vhich  occurs  during  a  quasi-steady  state 
period  for  each  cyclone,  is  sore  pronounced  in  cyclone  HI. 
This  cooling  within  each  budget  voluae  is  a  consequence  of 
the  development  and  poleward  movement  of  the  cyclones.  In 
addition,  although  the  troposphere  of  cyclone  N1  cools  more 
than  cyclone  S,  its  average  potential  teaperature  is  only 
slightly  lower  at  the  final  "data"  tine,  since  cyclone  HI 
was  initially  warmer  in  the  low  troposphere.  These 
potential  temperature  cross-sections  show  that  cyclone  N1 
initially  had  weaker  vertical  stability  than  cyclone  S. 
However,  by  the  end  of  the  study,  the  two  cyclones  have 
siailar  vertical  stability  characteristics. 


The  cooling  in  cyclone  N1  gives  the  appearance  of  a  cold 


ft 


« 
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core  strnctare.  This  however  is  not  the  case  3ince  the  cen¬ 
ter  of  each  cyclone  is  actually  located  in  the  wariest  air 
within  each  budget  volume  (see  ?ig.  11).  The  cooling  is 
likely  a  result  of  cold  advection  into  the  equatorward  half 
of  each  budget  voluie.  The  low- tropospheric  thermal  pattern 
can  be  seen  in  the  1000  mb  to  mid-layer  4  thickness  pattern 
for  hour  168  (see  Pig.  11).  The  thickness  pattern  shows  a 
warm  tongue  extending  from  the  east  through  the  center  of 
cyclone  Ml,  while  a  similar  warm  tongue  extends  east  to  west 
through  the  radius  nine  budget  volume  of  cyclone  S.  The 
thickness  pattern  also  reveals  that  the  cold  air  occupies 
more  of  the  budget  volume  of  cyclone  N1  than  cyclone  S,  thus 
resulting  in  the  slightly  lower  average  temeperatures  for 
cyclone  HI. 

Although  the  model  atmosphere  has  no  specified  tropo- 
pause.  Figs.  28  and  29  clearly  show  the  existence  of  a 
strong  tropopause  and  a  stratosphere  during  the  72  hours  of 
this  simulation  study.  Level  2  (approximately  200  mb) 
effectively  functions  similar  to  a  tropopause  layer  with 
slight  warning  and  strong  vertical  stability  above  200  mb. 
The  potential  temperature  evolutions  are  characterized  by 
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slight  fluctuations  belov  level  2.  These  fluctuations 
provide  a  hint  of  a  cyclic  evolution  pattern  during  the 
developaent  of  these  storms,  although  the  period  of  the 
cycle  is  not  evident. 

C.  HASS  TENDENCY 

There  are  no  sources  or  sinks  of  mass  within  the  atmos¬ 
phere,  however  it  is  possible  for  a  limited  volume  of  atmos¬ 
phere  to  gain  or  lose  mass.  The  time  rate  of  change  of  mass 
within  a  storm- following  budget  volume  provides  a  quantita¬ 
tive  measure  of  cyclogenesis.  Within  the  mass  budget,  the 
total  aass  within  each  specified  budget  volume  is  computed 
for  each  "data"  time.  From  these  values,  the  mass  tendency 
between  "data"  times  is  determined,  and  subsequently 
converted  to  a  pressure  tendency. 

The  aass  tendency  at  3 -hour  intervals  for  both  cyclones 
is  presented  in  Fig.  30.  A  cyclic  pattern  is  clearly 
indicated.  These  time  profiles  also  indicate  that  the 
developaent  of  the  model  cyclones  occurs  in  two  phases--an 
initial  explosive  cyclogenesis  period  (hours  96  to  126)  and 
a  quasi-steady  state  period  (hours  126  to  168).  A  12-hour 
oscillation  is  clearly  evident  throughout  these  two  phases. 
During  the  explosive  cyclogenesis  phase,  the  cycle  in  the 


radios  six  budget  consists  of  a  Bazinas  3-hour  pressure  fall 
followed  by  a  9-hour  period  of  relatively  snaller  surface 
pressure  change.  The  cyclic  pattern  during  the  quasi-steady 
state  phase  is  characterized  by  alternating  periods  of 
deepening  and  filling.  The  12-hour  cycle  can  be  seen  in  the 
pressure  tendency  patterns  for  both  the  inner  and  outer 
radii.  The  pressure  tendency  trace  for  each  budget  volume 
shows  that  thecyclone's  area-averaged  surface  pressure  drops 
in  discrete  intervals  like  a  step  function  rather  than  show¬ 
ing  continuous  deepening.  These  oscillations  could  be 
resoved  by  using  an  highpass  filter.  Piltering  would  also 
serve  to  emphasize  the  deepening  during  the  first  24  hours 
and  the  filling  which  occurs  later. 

The  pressure  tendency  oscillations  are  not  as  pronounced 
in  the  graphs  of  the  central  surface  pressures  of  these 
cyclones  (Fig-  31).  This  is  especially  true  during  the 
explosive  cyclogenesis  phase,  when  a  generally  steady 
decrease  is  seen  in  the  surface  pressure.  However,  the  rate 
of  change  nearly  corresponds  to  the  pressure  tendency  fluc¬ 
tuations  in  Fig.  30.  During  the.  second  phase  of  cyclone 
developient,  oscillations  are  seen  in  the  central  surface 
pressure,  however  they  are  not  as  large  as  those  the 


pressure  tendency 


The  cause  of  these  mass  fluctuations  is  not  known.  One 
possibility  is  that  the  cycle  is  the  reflection  of  an  exter¬ 
nal  gravity  wave  propagating  within  the  aodel.  An  external 
wave  would  have  a  siailar  effect  on  the  surface  pressures  of 
neighboring  anticyclones.  For  this  reason  a  study  of  the 
central  surface  pressures  of  aodel-generated  anticyclones 
was  completed. 

This  preliminary  study  of  five  high3  (three  in  the  HH 
and  two  in  the  SH)  for  periods  of  48  to  72  hours  (the  length 
of  the  period  was  dependent  on  whether  the  anticyclone  was 
an  analyzed  feature  at  a  particular  "data"  time)  strongly 
indicates  a  cyclic  pattern  of  surface  pressure  fluctuations 
(approximately  two  millibars)  in  the  model  anticyclones  as 
well.  The  pressure  oscillations  of  the  cyclones  and  anticy¬ 
clones  suggest  a  propagating  external  wave  within  the  model. 

The  pressure  tendency  traces  show  larger  fluctuations 
for  radius  nine  than  radius  six.  This  indicates  that  the 
pressure  falls/rises  for  the  area  between  radii  six  and  nine 
are  larger  than  those  for  the  area  within  the  radius  six 
budget  circle.  An  examination  of  other  budget  radii  showed 
that  the  magnitude  of  the  pressure  tendency  fluctuations 


increases  for  larger  radii 


This  also  suggests  a 
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large-scale  (perhaps  hemispheric  or  global)  sloshing. 
Because  the  scale  of  this  feature  is  auch  larger,  and  it 
does  not  interact  with  the  synoptic-scale  features,  the 
slosh  is  independent  of  cyclone  development.  A  aore 
detailed  investigation  of  the  aodel  surface  pressure  field 
is  required  to  isolate  this  phenomena  and  determine  if  it 
adversely  affects  the  numerical  simulations. 

D.  LATERAL  HASS  TRANSPORT 

One  goal  of  this  thesis  is  to  isolate  the  lateral  and 

vertical  branches  of  the  mass  circulation  of  the  model 

cyclones.  The  lateral  mass  transport  which  is  the  mass  flux 

through  a  specified  budget  volume  boundary.  The  lateral 

mass  transport  time  sections  (Pigs.  32  and  33)  graphically 

• 

describe  a  classical  cyclone  structure  of  low-level  converg¬ 
ence  (surface  to  approximately  level  4:  surface  to  600  mb) 
and  upper-level  divergence  (approximately  level  4  to 
mid-layer  2:  600  to  150  mb).  A  similar  pattern  is  noted  at 

radius  nine,  however  the  level  of  nondivergence  (LND)  is 
slightly  higher  (approximately  lid-layer  4:  500  mb)  .  The 

lateral  mass  transport  time  sections  indicate  that  the  mass 
circulation  does  not  extend  into  the  top  two  model  layers. 
This  further  substantiates  tha  existance  of  tropopause  and 
stratosphere  behavior  within  the  model. 
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The  aaxiaua  inflow  at  the  inner  and  outer  radii  of  each 


cyclone  occurs  in  the  lowest  aodel  layer  (centered  at 
approximately  900  ab) ,  while  th9  aaziwua  outflow  is  centered 
near  aid-layer  3  ( approximately  300  ab)  at  the  tropopause. 
In  addition,  the  aass  circulation  at  radius  nine  of  each 
cyclone  is  aore  intense  during  the  quasi-steady  state  phase 
than  during  the  rapid  cyclogenasis  phase.  This  is  a  further 
indication  that  the  the  second  phase  of  developaent  is  actu¬ 
ally  an  outward  growth  period  aore  than  a  deepening  period 
for  these  cyclones. 

i  cyclic  pattern  is  evident  at  the  levels  of  aaxiaua 
inflow  and  outflow  at  both  the  inner  and  outer  radii  of  each 
cyclone.  Looking  first  at  cyclone  HI  (Pig.  32),  the  aaxiaua 
outflow  at  radius  six  occcurs  at  hour  117  with  secondary 
aaxiaa  centered  at  hours  140,  150,  and  168  while  nearly 
equal  inflow  aaxiaa  are  seen  at  approximately  these  tiaes. 
The  aaxiaa  at  hour  117  occur  during  the  period  of  intense 
surface  developaent  (as  indicated  by  changes  in  the  central 
surface  pressure.  Pig.  31)  ,  while  the  maxima  at  later  tiaes 
occur  during  the  quasi-steady  state  phase  of  cyclone  Ml. 
The  lateral  aass  transports  at  radius;  nine  of  cyclone  Ml 
are  also  characterized  by  oscillations  in  the  basic 


inflow-outflow  pattern.  The  cyclic  variation  in  the 
lateral  aass  transport  are  coapled  with  the  oscillations 
in  the  pressure  field. 


ft  draeatic  increase  at  6°  and  a  decrease  at  9°  in  the 
height  of  the  level  of  LND  is  also  seen  during  the  quasi¬ 
steady  state  phase  of  cyclone  HI .  The  LND  height  change  for 
radius  nine  (approziaately  hour  13  5)  is  related  to  the 
intensification  of  the  aass  circulation  during  the  second 
stage  of  development  and  the  growth  of  the  cyclone  to  the  9° 
budget  volume  boundary. 

The  lateral  transport  tiae  sections  for  cyclone  S  also 
exhibit  the  classical  aass  circulation  pattern.  However, 
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the  cyclone  S  lateral  transport  pattern  is  noticeably  dif¬ 
ferent  from  that  of  cyclone  HI.  Specifically,  the  cyclone  S 
lateral  eass  transports  at  radios  3ix  are  more  intense  than 
those  of  cyclone  HI.  This  is  a  reflection  of  the  9reater 
intensity  of  the  spring  hemisphere  cyclone.  In  addition, 
the  height  of  the  LND  at  radios  six  of  cyclone  S  gradually 
lovers  fros  about  500  ab  to  near  700  ab,  and  is  generally 
lover  than  its  NH  counterpart  (which  remains  near  600  ab 
throughout  this  study)  after  hoar  132. 

The  radius  nine  tiae  sections  for  these  tvo  cyclones 
also  display  similarities  and  differences.  As  in 
cyclone  Hi,  there  is  a  simultaneous  decrease  in  the  height 
of  the  LND  and  the  intensification  of  the  lateral  aass 
transports.  However,  this  occurs  earlier  in  cyclone  S  (hour 
120  for  cyclone  S  and  hour  135  for  cyclone  N1) ,  and  the 
decrease  in  the  LND  height  is  aore  dramatic  (from  500  ab  to 
700  ab  for  cyclone  S  and  about  475  ab  to  600  ab  for 
cyclone  HI).  The  sharp  decrease  in  the  height  of  the  radius 
nine  LND  in  cyclone  S  occurs  sin ultaneously  with  the  outward 
growth  of  the  cyclone  to  the  outer  budget  volume  boundary. 
As  with  radius  six,  the  radius  nine  lateral  transports  are 
aore  intense  in  cyclone  S. 


Of  special  note  is  the  abrupt  weakening  of  the  cyclone  S 
aass  circulation  during  the  final  21  hours  of  this  study. 
This  weakening  is  most  noticeable  and  occurs  earlier  at 
radius  six.  The  decrease  in  the  aass  transport  intensity 
can  be  seen  in  both  the  inflow  and  outflow  branches,  however 
it  is  more  pronounced  in  the  outflow  branch.  This  weakening 
of  the  lateral  aass  transports  is  indicative  of  the  decay  of 
cyclone  S  and  is  reflected  in  a  corresponding  increase  in 
the  central  surface  pressure  of  cyclone  S.  In  contrast,  the 
lateral  aass  transports  at  the  cyclone  N1  budget  radii  con¬ 
tinues  to  be  strong  at  the  last  "data"  tiae  when  cyclone  N1 
reaches  its  lowest  central  surface  pressure  (984  ab) . 

E.  VERTICIL  HASS  TRANSPORT 

The  vertical  branch  is  an  integral  part  of  a  cyclone's 
aass  circulation.  Since  the  vertical  aass  transports  are 
the  area-averaged  vertical  velocity  structures  within  the 
cyclone  (in  both  isobaric  and  signa  coordinates) ,  tiae  sec¬ 
tions  of  the  aodel  vertical  velocities  (Pigs.  34  and  35)  are 
used  for  the  following  discussion  of  the  vertical  trans¬ 
ports.  This  discussion  points  out  siailarities  and  differ¬ 
ences  in  the  patterns  of  thess  two  stores.  In  general, 
upward  notion  is  seen  below  level  2  (approximately  200  mb) 
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at  the  inner  and  outer  radii.  Alternating  period,  of  weak 
upward  and  downward  aotion  ace  seen  above  level  2  in  the 
aodel  stratosphere.  The  upward  aotion  is  quite  intense  at 
radius  six  of  each  cyclone.  The  level  of  aaxiaua  upward 
vertical  notion  coincides  with  the  LMD  indicated  in  the  lat¬ 
eral  aass  transport  tine  sections  (Figs.  32  and  33). 
periods  of  aaxiaua  upward  vertical  notion  coincied  (as  they 
aust)  with  the  periods  of  aost  intense  lateral  aass  trans¬ 
port.  transport.  The  upward  vertical  notion  at  radius  nine 
is  weaker  than  at  radius  six.  However,  the  intense  upward 
vertical  aotion  occurs  only  within  the  troposphere. 

Considering  first  the  six  degree  budget  volune,  two  dis¬ 
tinct  periods  of  aaxiaua  upward  vertical  notion  are  evident 
in  each  cyclone.  For  cyclone  HI ,  the  first  period  (hour  107 
to  123)  occurs  during  the  explosive  cyclogenesis  phase, 
while  the  second  period  (after  hour  147)  occurs  during  the 
qi.'&3i-steady  state  phase.  The  subsequent  radius  six  aaxinun 
occurs  siaultaneously  with  a  siailar  aaxiaua  at  radius  nine, 
and  results  froa  the  outward  growth  of  the  cyclone  to  the 
outer  budget  voluae  boundary. 

Cyclone  S  also  experiences  two  periods  of  intense  upward 
vertical  aotion  (centered  at  hours  117  and  138).  As  in 


cyclone  Ml,  the  first  aaxiaua  occurs  during  the  explosive 
cyclogenesis  phase  vhile  the  second  aaxiaua  and  its  radius 
nine  counterpart  occur  during  the  period  of  outward  lateral 
growth.  Hear  the  beginning  of  the  integration,  a  third 
period  of  intense  upward  vertical  motion  is  evident  for  each 
budget  radius.  However,  its  significance  is  unclear  due  to 
the  lack  of  the  proceeding  model  data.  In  addition,  the 
decay  of  the  cyclone  is  indicated  by  the  rapid  decrease  in 
upward  vertical  motion  after  hour  159. 

in  interesting  feature  of  the  cyclone  S  radius  nine  tine 
section  is  the  downward  vertical  motion  indicated  in  the 
lower  two  layers  between  hours  101  and  111.  This  requires  a 
downward  notion  within  the  6°-9°  area  which  more  than  off¬ 
sets  the  upward  motion  within  radius  six.  This  indicates 
strong  subsidence  around  the  cyclone  compensating  for  the 
intense  upward  vertical  motion  associated  with  the  explosive 
cyclogenesis  phase  when  cyclone  S  is  still  within  the  6° 
budget  volume.  This  is  a  possible  reflection  of  the  cyclic 
evolution  of  the  model  cyclones. 

Time  traces  of  area- averaged  vertical  velocities  for 
level  4  of  cyclones  HI  and  S  are  presented  in  Pig.  36.  The 
two  major  periods  of  strong  upward  vertical  motion  which 
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correspond  to  the  two-phase  cyclone  development  are  clearly 
depicted.  The  outward  growth  phase  of  each  cyclone  is  viv¬ 
idly  displayed  in  these  traces.  A  draeatic  increase  in  the 
radius  nine  area-averaged  vertical  velocities  is  seen  as  the 
lateral  growth  of  these  cyclones  exceeds  the  9°  budget 
voluae  boundary. 

The  radius  six  trace  for  cyclone  If  is  characterized  by 
two  distinct  aaxiaa  at  hours  115  and  168.  The  radius  nine 
trace  essentially  has  a  single  saxiaua  which  occurs  during 
the  second  phase  of  cyclone  developaent.  The  extreae  aax- 
iaua  at  hour  115  is  approxiaately  in  the  aiddle  of  the  rapid 
cyclogenesis  phase  when  the  central  surface  pressure 
decreases  17  ab  between  hours  13  8  and  123  (see  Pig.  31). 

A  surprising  feature  of  the  vertical  velocities  occurs 
between  hours  141  and  150  when  the  aagnitude  within  radius 
nine  actually  exceeds  the  corresponding  value  for  radius 
six.  This  inplies  that  the  vertical  velocities  within  the 
go. go  region  exceed  those  within  the  radius  six  budget 
circle. 

Although  the  cyclic  pattern  seen  in  cyclone  S  is  siailar 
to  that  of  cyclone  9,  there  are  significant  differences.  As 
in  cyclone  N,  two  aaxiaa  are  seen  at  radius  six 
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(hoars  118  and  136),  however  the  tiee  interval  between  these 
two  occurrences  is  each  shorter  than  in  cyclone  H.  The 
first  aaxiaue  occurs  daring  the  first  phase  of  cyclone 
developeent  when  the  central  surface  pressure  decreases 
22  eb  between  hours  99  and  126.  The  second  aaxiaua,  which 
occurs  during  the  quasi-steady  state  phase  of  cyclone  S, 
coincides  with  the  single  aaxiaua  at  radius  nine.  The  aaxi- 
aue  at  radius  nine  of  cyclone  S  is  larger  than  its  cyclone 
H  counterpart.  The  decay  of  cyclone  S  is  clearly  indicated 
by  the  sharp  decrease  in  the  level  4  vertical  velocities  at 
both  six  and  nine  degrees. 

?.  COBP A BIS OH  OF  MODEL  AND  BUDGET  VERTICAL  VELOCITIES 

One  objective  of  this  thesis  is  a  coaparison  of  the 
aodel  vertical  velocities  with  those  coaputed  in  the  aass 
budget.  The  vertical  divergence  of  the  vertical  aass  trans¬ 
port  (budget  residual)  and  the  aass  tendency  are  coaputed 
within  the  aass  budget  (see  Tables  IV  and  V).  A  vertical 
integration  of  these  values  (fro a  an  upper  boundary  value  of 
a  0  at  »  0)  describes  the  vertical  aass  transport,  which 
is  then  converted  to  an  area- averaged  vertical  velocity. 
Thus  the  budget  vertical  velocities  are  essentially  kine- 
aatic  oaega  values  for  the  specified  budget  voluae. 


M  10t  10*  1M  tao  131  132  1M  1 
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Figure  36.  flodel  vertical  velocity 
(-//bar  sec-i)  a.  ana 


This  approach  has  been  used  in  atsospheric  studies  to 
infer  vertical  velocities  and  vertical  transports  of  nass 
and  other  properties.  In  this  study,  vertical  velocities 
are  available  as  part  of  the  nodel  "data"  base.  This 
affords  the  unique  opportunity  to  test  the  accuracy  of  the 
budget-derived  vertical  velocities. 

The  basic  patterns  in  the  nodel  vertical  velocity  tine 
sections  (Pigs.  34  and  35)  are  contained  in  the  budget  ver- 
tical  velocity  tine  sections  (Pigs.  37  and  38) .  The  budget 
vertical  velocities  for  radii  3ix  and  nine  of  each  cyclone 
are  aost  intense  belov  the  tropopause  (approximately  mid¬ 
layer  3)  .  Very  weak  vertical  notion  is  seen  in  the  nodal 
stratosphere,  At  radius  six  of  each  cyclone,  two  maxina  are 
seen  in  the  budget  vertical  velocities  tine  sections,  with 
the  earlier  aaxinua  occurring  during  the  first  phase  of 
cyclone  development.  is  in  the  nodel  vertical  notion  pat¬ 
tern,  the  second  aaxinua  in  the  budget  vertical  velocity 
tiae  section  occurs  in  conjunction  with  a  radius  nine  maxi¬ 
nun.  This  again  indicates  the  lateral  growth  of  the 
cyclones  beyond  the  9°  budget  volume  boundary.  For  each 
budget  radius  of  cyclone  S  there  is  a  significant  decrease 
in  the  budget  vertical  velocities  after  hour  159  when 


cyclone  S  begins  to  fill.  This  feature  is  also  seen  in  the 
■odel  vertical  velocities. 

Although  the  patterns  in  the  eodel  and  budget  vertical 
velocity  tiae  sections  are  siailar,  significant  differences 
vere  noted  in  the  aagnitudes  of  the  budget  vertical  veloci¬ 
ties.  Tiaesections  of  the  differences  between  these  two 
values  (budget  value  ainus  nodal  value)  are  presented  in 
Pigs.  39  and  40.  With  the  exception  of  the  stratosphere, 
the  budget  vertical  velocities  were  larger  than  those  of  the 
aodel.  It  appears  that  the  largest  departures  occur  where 
the  vertical  velocities  (both  aodel  and  budget)  are  largest. 

Since  the  aodel  has  an  internally  consistent  "data"  base 
there  are  two  possible  sources  for  the  errors  in  the  budget- 
derived  vertical  velocities.  Errors  can  result  froa  the 
interpolation  of  the  winds  froa  the  rectangular  grid  of  the 
aodel  to  the  circular  grid  of  the  aass  budget.  Other  errors 
arise  froa  the  vertical  distribution  of  the  inbalances 
between  the  integrated  divergence  and  the  aass  tendency, 
which  aust  balance  since  there  are  no  sources  or  sinks  of 
aass.  The  coaparison  of  the  budget  and  aodel  vertical 
velocities  indicates  that  the  lateral  aass  transports  are 
too  large  after  the  interpolation  and  balancing  proceedures 


are  completed  within  the  aass  budget.  Given  the  consistent 
"data"  of  the  eodel,  these  errors  are  generally  larger  than 
one  eight  anticipate.  1  detailed  exaeination  of  the  inter¬ 
polation  scheee  and  the  vertical  distribution  of  the  budget 
residuals  should  be  conducted  in  an  effort  to  isolate  the 


source  of  these  errors 
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71.  CONCLUSIONS 


In  this  thesis  two  numerically  simulated  oceanic  extra- 
tropical  cyclones  were  studied.  Initially  a  synoptic  study 
of  these  cyclones,  including  a  comparison  with  observed 
cyclones,  was  completed.  These  cyclones  were  four  to  be 
characterized  by  both  realistic  and  unusual  synop  *  fea¬ 
tures.  In  a  general  sense,  the  realistic  featu  were 
associated  with  the  lower  troposphere  while  the  unusual,  fea¬ 
tures  were  associated  with  the  middle  and  upper  troposphere 
and  the  vertical  coupling  of  the  model  layers. 

These  cyclones  also  vere  studied  through  the  application 
of  guasi-Lagrangian  diagnostics.  This  diagnostic  study  of 
the  mass  circulation  of  model  cyclones  showed  a  concentra¬ 
tion  of  inward  lateral  transport  (convergence)  in  the  lowest 
model  layer  and  outward  lateral  transport  (divergence)  in  a 
layer  centered  near  300  mb.  Although  the  model  has  no  spec¬ 
ified  tropopause,  the  basic  patterns  of  the  lateral  and 
vertical  mass  transports  and  the  thermal  structure  revealed 
the  existence  of  a  tropopause  near  mid-layer  3  of  the  model. 

Traces  of  the  mass  tendency  within  6°  and  9°  budget 
volumes  revealed  a  12-hour  cyclic  pattern  imbedded  in  a 
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two- phase  cyclogenesis  pattern.  The  first  phase  of  cyclone 
development  was  a  period  of  sc  plosive  deepening  while  the 
second  phase  was  a  quasi-steady  state  period.  The  12-hoar 
cycle  was  aost  easily  seen  in  the  pressure  tendency  trace. 
Daring  the  first  phase  of  cyclone  development,  the  cycle 
consisted  of  a  maximum  3-hour  pressure  fall  followed  by  a 
9- hour  period  of  little  surface  pressure  change.  During  the 
second  phase,  the  cycle  was  characterized  by  alternating 
pericds  of  pressure  falls  and  rises.  It  was  also  noted  that 
the  fluctuations  in  the  pressure  tendency  increased  for  suc¬ 
cessively  larger  budget  radii.  The  two-phase  cyclone 
development  was  seen  in  the  traces  of  the  central  surface 
pressures.  However,  the  12-hour  cyclic  pattern  was  not  eas¬ 
ily  detected  in  these  traces.  Phis  12-hour  cycle  produces  a 
step  function  decrease  in  the  area-averaged  surface  pressure 
of  these  cyclones  rather  than  a  steady  deepening.  This  step 
function  evolution  is  not  clearly  indicated  by  the  central 
surface  pressure  traces  (see  Pig.  31).  &  possible  cause  of 
this  cycle  is  an  external  wave  propagating  within  the  model. 
However,  further  studies  are  needed  to  isolate  the  source  of 
these  oscillations. 
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The  lateral  aass  transport  and  the  vertical  aass  trans¬ 
port  (converted  to  a  vertical  velocity)  clearly  described 


the  aass  circulation  of  these  cyclones.  Tiae  sections  of 
these  terns  also  indicated  two-phase  cyclone  develpaent. 
Two  aaxiaa  were  seen  at  radius  six  of  each  cyclone— one 
during  each  phase  of  developaent .  The  one  aaximua  seen  at 
radius  nine  occurred  in  conjunction  with  the  second  aaxiaun 
at  radius  six  and  indicated  the  outward  growth  of  the 
cyclones  beyond  the  outer  budget  volune  boundary.  Haxiaa  in 
the  lateral  aass  transport  and  vertical  velocity  tiae  sec¬ 
tions  generally  occurred  siaulta neously.  In  addition,  the 
onset  of  the  decay  of  cyclone  5  was  indicated  by  the  rapid 
decrease  in  the  intensity  of  the  lateral  aass  tendency  and 
vertical  aotion  near  the  end  of  the  study. 

Area-averaged  vertical  velocities  were  *  computed  within 
the  aass  budget  and  were  coapared  to  area-averaged  model 
vertical  velocities.  Although  the  tiae  sections  of  these 
terns  have  siailar  basic  patterns,  significant  differences 
were  found  in  their  magnitudes.  It  appeared  that  the  depar¬ 
tures  were  largest  where  the  vertical  velocity  (both  model 
and  budget)  were  the  largest.  The  two  possible  sources  of 


these  errors  are  interpolation  of  the  wind  from  the  model 
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grid  to  the  budget  voluae  grid,  and  iabalances  in  the 
vertical  distribution  o£  the  aa3s  budget  residuals.  Further 
studies  are  required  to  isolate  the  relative  magnitudes  of 


these  two  errors 
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The  results  of  the  synoptic  study  and  the  sass  budget  analy¬ 
sis  suggest  that  further  diagnostic  studies  of  the  basic 
properties  of  these  cyclones  hre  needed  to  fully  understand 
their  dynasic  structure.  In  addition,  further  studies  are 
needed  to  determine  the  source  of  the  errors  in  the  budget 
derived  vertical  velocities,  with  this  in  eind,  the  follow¬ 


ing  recossendations  are  aade: 


•  Exaaine  the  budget  technique  for  coaputing  area- 
averaged  vertical  velocities  and  the  interpolation 
scheae  eaployed  within  the  aass  budget  to  isolate  the 
sources  of  the  errors  in  the  budged  derived  vertical 
velocities. 


•  Per fora  budget  studies  on  the  aodel  cyclones  during 
earlier  and  later  tiae  periods  than  covered  in  this 
study. 


•  Per  fora  budget  studies  on  aodel  generated  anticyclones 
to  deteraine  if  they  are  characterized  by  the  l 2-hour 
cyclic  pattern  of  the  aodel  cyclones. 


•  Study  the  cyclic  oscillations  by  differencing 
successive  pressure  fields. 


•  Perfora  budget  studies  on  cyclones  and  anticyclones 
generated  in  the  Sandgatne  (1981)  adiabatic 
experiments. 


•  Perfora  budget  studies  of  other  Droperties  on  these 
aodel  cyclones  and  on  aodel  anticyclones. 


APPENDIX  & 


COMMONLY  USED  ACRONYMS  AND  ABBHE YIATIONS 

NOSAPS  Navy  Operational  Ataoapheric  Prediction  Model 
UCLA  University  of  California,  Los  Angeles 
NEW  Pacxiity*^CO&aent&^  Prediction  Research 

SHE  east-northeast 

ESE  east-southeast 

NNE  north- northeast 

SSR  south- southwest 

LND  *  level  of  nondivergence 
AVA  anticy clonic  vorticity  advection 

CVA  cyclonic  vorticity  advection 

PBL  planetary  boundary  layer 

sec*1  per  second 
■  b  sillibar 

°C  degree  Celcius 


° R  degree  Kelvin 


SH  northern  hemispheres 

SB  southern  hemisphere 

ns  nautical  miles 

B  mass 

S  north 

S  south 

a  meters 
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APPBWDIX  B 

QOASI-LAGHAMGIA E  BUDGET  FRAHEIDBR  AHD 
GEMEBALIZED  BUDGET  EQUATIONS 

The  purpose  of  this  appendix  is  to  present  the  foraula- 
tion  of  the  store  budget  voluae  coordinate  systee  and  the 
generalized  budget  equations  which  can  be  applied  to  any 
desired  ataospheric  property.  The  inforaation  contained 
herein  is  extracted  froa  Bash  (1978)  as  adopted  froa  Johnson 
and  Dovney  (1975  a). 

The  stora  budget  voluae  in  spherical  coordinates  is 
depicted  in  Fig.  41.  The  following  notation  is  relevant  to 

Aziauthal  coordinate 

Angular  radial  coordinate 

Position  vector  froa  earth's 
center 

position  vector  froa  stora  center 
to  any  point  in  budget  voluae 

Orthognal  unit  vectors  along 
aeridional.  aziauthal  and, 
vertical  directions  of  spherical 
coordinates  for  stora 

find  velocity  relative  to  earth 


this  figure: 
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» 


Horizontal  velocity  of  reference 
axis 


90  Horizontal  velocity  of  any  point 

~  on  the  volaee  boundary. 


Table  VI  contains  a  listing  of  the  generalized  budget 
equations.  The  following  is  an  explanation  of  the  notation 
used  in  these  equations: 


a  The  vertical  coordinate  systea 

V  Toluse 


T  Top  boundary 

B  Bo t toe  boundary 


Horizontal  wind  coaponent  relative 
to  budget  voluae  noraal  to  lateral 
bounds  ry 


(d<r/dt)  Vertical  velocity  divergence  into 
a  layer  within  the  voluie 


rain/? 

p* 

f 

b 


Radius  of  budget  voluae 
Surface  pressure 

Specific  budget  property  (">er  unit  aass) 
Vertical  velocity  (d<r/dt) 


9 


Acceleration  duo  to  gravity 
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TABLE  V 

GENERALIZED  QUASI-LAGRA NGI AN  BUDGET  EQUATIONS  IN  SIGHA 

COORDINATES 
(after  Nash ,  1978) 


The  budget  woluse  integral  is  defined  as 
2 T  r  B  2r  p*  2 

p  a  I  I  I  — — -  fr  ainBdcidBC-da)  • 

°B  0  0 


The  east  general  fora,  of  the  budget  equation  is 

+  VT(F)  +  S(F) 


£  “  »  <*> 


where  the  lateral  transport  tern  is 
r^r  2* 


LT(F) 


ff  r  p* 

J  I  -  f(D-W)  &r  sin  8do(-da)  |&B  , 

°B  J0  * 

the  vertical  transport  tern  is 

Si 


VT(F) 


ft  "B  f 

k  L  ~r  h <fo>  * 


ingdodB  (-do) , 


'B  0  '0 

and  the  source  (sink)  tern  is 

?t  2* 

SW  "J,  J  f  T^'! 

BoO 

The  wertical  transport  tern  is  zero 
when  the  lisits  of  integration  are 
o*l  and  o*0y  since  there  can  be 
no  transport  through  the  top  and 
bottoa  boundaries  of  the  budget  woluse* 
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